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Abstract 

In this paper, an extended Kalman filter (EKF) 
is used to estimate the parameters of a low-order 
model from aircraft transient response data. The low- 
order model is a state space model derived from the 
short-period approximation of the longitudinal air- 
craft dynamics. The model corresponds to the pitch 
rate to stick force transfer function currently used 
in flying qualities analysis. Because of the model 
chosen, handling qualities information is also ob- 
tained. The parameters are estimated from flight 
data as well as from a six-degrce-of-frecdom, non- 
linear simulation of the aircraft. These two esti- 
mates are then compared and the discrepancies noted. 
The low-order model is able to satisfactorily match 
both flight data and simulation data from a high- 
order computer simulation. The parameters obtained 
from the EKF analysis of flight data are compared 
to those obtained using frequency response analy- 
sis of the flight data. Time delays and damping ra- 
tios are compared and are in agreement. This tech- 
nique demonstrates the potential to determine, in near 
real time, the extent of differences between computer 
models and the actual aircraft. Precise knowledge 
of these differences can help to determine the flying 
qualities of a test aircraft and lead to more efficient 
envelope expansion. 


•Aerospace Engineer. Member AIAA. 

••Professor. Member AIAA. 

Copyright ©1990 by the American Institute of Aeronautics and 
Astronautics, Inc. No copyright is asserted in the United Slates 
under Title 17, U.S.Code. The U.S. Government has a royalty- free 
license to exercise all rights under the copyright claimed herein 
for Governmental purposes. All other rights are reserved by the 
copyright owner. 


Nomenclature 

d distance from angle-of-attack vane to aircraft 

center of gravity 

EKF extended Kalman filter 

F(x,t ) linearization of / about % at time t 

f(x,t ) nonlinear function dependent on 

variable x at time t 

g gravitational constant 

H linearized output matrix, defined in 

equation (22) 

h nonlinear output matrix 

I identity matrix 

K Kalman gain matrix 

L q derivative of a with respect to q 

L a derivative of a with respect to a 

Ls derivative of a with respect to 5 

LOES low-order equivalent system 

M q dimensional variation of pitching 

moment with pitch rate 

M a dimensional variation of pitching 

moment with angle of attack 

Ms dimensional variation of pitching 

moment with surface deflection 

n* normal acceleration 

P Riccati error covariance 

Q model noise variance (process noise) 

q pitch rate 

R measurement noise variance 


V 

velocity 

V 

measurement noise 

X 

state vector 

a*(+) 

estimate of x taking advantage 
measured data at k 

£*(-) 

estimate of x before measured 
data at k is available 


To provide more detailed information, an additional 
real-time technique, the subject of this paper, has been 
investigated. This technique uses an extended Kalman 
filter (EKF) to estimate a low-order equivalent sys- 
tem (LOES) model of the aircraft from flight-measured 
time histories. References 4-9 detail other work in 
this area. For other work in parameter estimation see 
Refs. 10 and 11. 


z output vector 

at angle of attack 

8 pseudo control surface deflection 

S p pilot pitch stick signal 

p correlation coefficient 

o standard deviation 

t first-order time lag coefficient 

Td pure time delay 

w process noise 

short-period natural frequency 

Subscripts 


/ 

k 

m 

0 


corrected 
flight measured 
at discrete time k 
measured parameter 
initial 


Superscript 

T matrix transpose operator 

estimate 


Introduction 

Many near-real-time techniques have been used in 
the envelope expansion of recent research aircraft like 
the X-29. 1 In the area of flight controls, two real-time 
techniques were developed. One is a frequency re- 
sponse technique 2 where the aircraft’s open-loop fre- 
quency response is determined from measured flight 
parameters. This frequency response is then compared 
with a predicted frequency response based on the cur- 
rent aircraft simulation. The second technique is a time 
history overplot comparison, 3 conducted in real time. 
Overplots are obtained of the aircraft response with 
time histories generated from a previously determined 
linear model of the aircraft, at that flight condition, 
driven by the pilot’s input signal as it is telemetered 
to the ground. 


A Kalman filter is an optimal state estimator of 
a linear system. The estimate is adjusted based 
on each measurement of the system. An extended 
Kalman filter extends linear Kalman filter theory to 
nonlinear systems. 

The LOES model is based on the longitudinal short- 
period approximation in the MIL-STD-1797 for fly- 
ing qualities. 12 This technique is able to extract fly- 
ing qualities parameters in near real time, which may 
quantify discrepancies between the aircraft and the air- 
craft simulation. 

Before the flight, the aircraft’s six-degree-of- 
freedom, nonlinear simulation is analyzed at a speci- 
fied flight condition, and a LOES is estimated. During 
flight test, parameters of the low-order model are es- 
timated and these flight-data based estimates are com- 
pared with the simulation based estimates in real time. 
The results of the flight estimation are intended to help 
diagnose discrepancies between the aircraft response 
and the predicted response. This comparison could be 
done before the aircraft reaches its next test point, giv- 
ing the flight controller time to decide whether the dis- 
crepancy is serious enough to curtail further testing. 

The aircraft used in developing this technique was 
the X-29 (Fig. 1). The X-29 airframe has a negative 
static margin of approximately 35 percent. This in- 
stability requires a high degree of control augmenta- 
tion. The control system is entirely fly-by-wire, con- 
sisting of a triply redundant primary channel operating 
at 80 Hz. The longitudinal control laws use pitch rate 
and normal acceleration feedbacks for primary stabil- 
ity and control. Of particular importance to the EKF 
algorithm is the airplane’s three surface pitch control, 
using canards, symmetric flaps, and strake flaps. Data 
from the airplane is telemetered to the control room at 
40 samples/sec (40 Hz). 

The estimation was made of the closed-loop system 
rather than the open-loop system. Because of the high 
degree of control augmentation on the X-29, the re- 
searchers thought that estimating the closed-loop sys- 
tem offered the best chance of success. The control 
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system was designed to make the closed-loop aircraft 
behave much like the low-order model. Detailed post- 
flight analysis would provide greater insight into any 
discrepancy by separating the effects of the control 
system from the aerodynamics. 

If a major discrepancy between the simulation and 
actual airplane is present, this algorithm may be able to 
determine what the discrepancy is more quickly than 
currently available technology can. 

Not only general information about the aircraft per- 
formance, compared to its predicted performance, is 
available, but handling qualities data is obtained as 
well. First, an equivalent time delay is available. This 
value can be monitored to identify potential stability 
problems. The time delay computed in near real time 
would be used to identify large discrepancies from the 
anticipated values. 

In addition, plots of normal acceleration divided by 
angle of attack (n* / a) could be obtained for both the 
simulation and real-time estimations. This would al- 
low researchers to determine if the flight values are 
moving towards an unacceptable region as defined by 
Ref. 12. 

With the information provided by this technique in 
near real time, more intelligent decisions about the on- 
going flight test could be made, possibly saving time 
and money. The results obtained by this technique 
could also increase the efficiency and safety of the 
flight test program. 

The potential for reducing postflight data process- 
ing is also great. This technique could determine the 
area responsible for any discrepancies, and postflight 
analysis could then identify the cause. 

Techniques and Procedures 

Modeling 

Only the transient longitudinal dynamics of the air- 
craft in l-p level flight are studied. The basic model is 
the longitudinal short-period approximation 

a = L a cn+ L q q + L&6 (1) 

q = M a ct + M q q + Af$6 (2) 


tion to be obtained. Aircraft angle of attack, pitch rate, 
and normal acceleration make up the systems observa- 
tion equations. 

Because of the aircraft’s open-loop instability, the 
closed-loop dynamics were estimated rather than the 
open-loop dynamics. If the open-loop equations were 
used, then any small measurement errors in the surface 
positions or the numerical differentiation would cause 
the model to diverge from the measurements. 

The aircraft uses three surfaces for pitch control. 
This configuration made defining a stick force to air- 
craft response transfer function difficult. Since there 
was no single pitch generating surface, the effects of 
the three surfaces were collected into a single pseudo 
control surface. A third state equation is added to the 
model to represent the pseudo control 

6(t) = —6(0 + — 6 p (t - r d ) (3) 

T T 

where 6 is the pseudo control surface deflection and 
6 P is the pilot pitch stick signal. It was reasoned that 
this equation should have the appearance of an actu- 
ator model. The coefficient of the model is estimated 
by the filter. Equation (3) contains a first-order lag of 
the pilot’s stick signal. The first-order lag accounts for 
only part of the system time delay. The rest is modeled 
through a pure time delay. 13 

To determine the pure time delay, several a priori 
estimates are made, and several filters are run in par- 
allel. To determine which of these filters is operating 
with the proper amount of pure time delay, the like- 
lihood function is used. 14 The proper amount of pure 
lime delay is the sum of each filter’s time delay multi- 
plied by its likelihood. However, in all cases examined 
one delay always resulted in a likelihood much higher 
than the others. For this reason, and to save compu- 
tational complexity, the time delay associated with the 
largest likelihood is used as the correct pure time delay. 
Figure 2 shows that the pure time delay is implemented 
outside of the EKF. 

A detailed investigation of identifiability was not 
conducted. However, the fact that one delay resulted in 
a likelihood far larger than the others suggested iden- 
tifiabilty is not a problem. 


where L„, L q , and Ls are the derivatives of a with Since the aircraft cannot obtain perfect trim in flight, 

respect to a, q, and 6, and M a , M q , and Mi are the estimated bias terms are added to the state equations, 

derivatives of q with respect to a, q, and 6. This ap- In the presence of trim errors these terms balance the 

proximation is equivalent to the short-period approxi- state equations and accurately determine the trim state, 

mation for the pitch rate to stick force transfer function Without the biases, the equations would integrate a 

given in Ref. 12, allowing handling qualities informa- small error over the course of the maneuver. 
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The equation for measured angle of attack needs to 
be corrected for the fact that angle of attack is mea- 
sured at the noseboom of the aircraft and not the cen- 
ter of gravity. This can be modeled to the first order by 
correcting for pitch rate 

a m = a-^q (4) 

where or TO is the measured angle of attack and d is the 
distance from the angle-of-attack vane to the aircraft 
center of gravity. 

The equivalent system model can now be defined 
and is given by equations (5-10) 


a = L a a + L q q + Ls& + Lo 

(5) 

q = M a oi + M q q + M&8 + Mo 

(6) 

8 = — 8 + — 6 P + <5o 

(7) 

r r 

d 

at m = a - y<! 

(8) 

►Cj 

3 

It 

(9) 

-V 

rh m - [L a a+ (Lq - l)g 

9 

+ Ls6 + f/ol + «zo 

GO) 


where the subscript 0 indicates initial state. 


The continuous system model is 

x(t) = /( x(t),0+ «>(t) (12) 

where w is a white noise process with spectral density 
Q(t) where Q is the model noise variance. 

The discrete measurement model is 

zk= h k ( x(t k )) + I/* ; k= 1,2, 3, ... (13) 

where the v k 's are independent Gaussian noise vectors 
with mean zero and covariance R k . It is also assumed 
that uj and v are independent. 

The initial conditions are 

i(0) (14) 

where x(0) is Gaussian, with mean xo and vari- 
ance Pq. 

The state estimate propagation between samples is 

iit) = /(*(*), o (i5) 

The error covariance propagation between sam- 
ples is 

P(t) = F(x(t),t)P(t) 

+ P(t)F J (x(t),t) + Q(t) (16) 


Estimation 

The estimation problem consists of estimating L a , 
L q , L s , M a , Mq, Ms, i, Lo, Mo, 8 0 , n* 0 of the lon- 
gitudinal short-period approximation model as a func- 
tion of time for the closed-loop aircraft response, in 
near real time. To do this using an EKF, these coeffi- 
cients must be made states of the filter. Therefore the 
state vector ( x) becomes 

X [ Oi t t] y & , L a y Kfa y L$ y Lf ft y Lj q y 

Mq, Lo, Mo, So, n* 0 ] T (11) 

The extra state equations required are L a = 0, M a = 
0, L s = 0. Ms = 0, L, = 0, Mq = 0, ( J) = 0, Lq = 0, 
Mo = 0, 6o = 0, and = 0. The short-period state 
equations (5-7), augmented by those previously listed 
and output equations (8-10), form the state and mea- 
surement equations for the EKF. 

The EKF equations used are for a continuous- 
discrete system and are listed in the following 
equations. 15 


where P(t ) is the Ricatti error covariance, and 
F( x(t) , 0 is defined in equation (20). 

The state estimate update is 

£*(+) = **(-) + K k [z k - h k (x k (-))} (17) 

where (+) is the estimate using information at time 
k, (— ) is the estimate before data from time k is 
available, and K is the Kalman gain matrix defined in 
equation (19). 

The error covariance update is 

P k (+) = [I-K k H k (x k (-))]Pk(-) (18) 

where H k (x k (-)) is defined in equation (21). 

The Kalman gain matrix is 

K k = P k (-)Hj(x k (-)) 

[H k (x(-))P k (-)Hj 
(x k (-)) + R k r x (19) 
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Definitions 


F(i(t),f) 


H k (i k (-)) 




(20) 

dx{t) 

*(0=a<‘) 

dhkWk)) 


(21) 

dx(t k ) 

*(<*)=**(-) 


Equation (12) for the system model is a continu- 
ous nonlinear equation describing the state propaga- 
tion with time, between samples. An approximation 
that F is constant between samples has been made. 
In this application, the state is propagated through the 
nonlinear equations using a fourth-order Runge-Kutta 
technique. The measurement model (equation (13)) 
is discrete. In equation (20), F is the linearized state 
matrix. This matrix must be recomputed for each up- 
date of the state vector since the elements of F are also 
states of the EKF and will change with each update of 
the filter. 

The first step in estimating the flight data is obtain- 
ing an estimate from the simulation. There are two 
primary objectives in estimating the low-order model 
from nonlinear simulation data. First, determine the 
coefficients of the low-order model (filter states 4-11) 
as accurately as possible. These will be the initial state 
estimates for the estimation from the flight data. In 
addition, the simulation estimates will be compared to 
the estimates obtained from the flight data to check for 
major discrepancies. 


signal. These bounds are then selected as ±2 a (two 
standard deviations) bounds about the mean. The 
value for Ri t i in this case is of *. This is possible be- 
cause the process noise was kept low. 

For a and q, the initial state variance ( P( 0) ) is the 
value for measurement noise on a and q. For the other 
slates a value for a is selected which equals one-fourth 
the estimated initial values of the states, thus Pi ( i = 
&? f = ( jXi'i ) 2 . For state variables with initial values 
close or equal to zero, bounds are established based on 
engineering judgement. 

With the identification algorithm initialized, the 
next step is to determine the pure time delay. Nine 
filters are run with fixed time delays ranging from 0 to 
200 msec, in increments of 25 msec. The likelihood 
function is used to determine which filter is operating 
with the correct amount of time delay. 

Several analysis techniques are used to monitor the 
filter and to aid in interpreting the results. The most 
important of these are the filter innovations. The inno- 
vations are the difference between the predicted output 
of the equivalent system and the measured response. 
The innovations are calculated by 

innovations - Zk- hk( x k( -)) ( 22) 

Results and Discussion 


Simulation 


The second objective is to determine the amount of 
process noise to use in estimating the data from flight. 
For the X-29A, earlier flight tests have shown that the 
nonlinear simulation accurately predicts the response 
of the airplane. As a result, the process noise weight- 
ing matrix is selected by adjusting the weighting ma- 
trix until both the EKF and LOES indicated that good 
results had been obtained. Ovcrplots of the LOES sys- 
tem with the original data are used to evaluate the esti- 
mates qualitatively. Innovations, time histories of the 
estimates, correlations, and state variances from the 
EKF are used in the evaluation. A tradeoff is then 
made between the EKF and LOES results to obtain the 
final estimates and weightings. These weightings are 
then used to estimate the parameters during flight. In 
the cases examined, the process noise weightings were 
kept quite low. 16 

To determine values for the measurement noise 
(R), sample time histories are obtained. From these, 
bounds are drawn defining a noise envelope for that 


A simulation time history is generated using the 
six-degree-of-freedom, nonlinear simulation and flight 
data inputs at this condition. Some of the flight data 
analyzed is presented in Fig. 3. The actual time inter- 
val of the identification was varied to illustrate vari- 
ous points. 

The filter is initialized as described earlier. The mea- 
surement noise is determined from the flight data and 
then divided by 100 (to represent the root mean 
square error) since the simulation would produce a per- 
fect measurement signal. A low-order equivalent sys- 
tem is then estimated from this data. The results of this 
estimation are presented in Table 1. 

Large correlations between the state 6 and the co- 
efficients Mo and <5o are found. Some correlation is 
also indicated between the coefficients Mq and M&. 
This is expected because the control system has the ap- 
pearance of a pitch rate command system. Although 
this correlation is not desirable, it is not felt to be 
a problem. 
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These results are obtained for very low levels of pro- 
cess noise. Large values of process noise often result 
in misleading estimates. The values for final variance 
are much smaller than the initial values, which shows 
that the EKF has obtained a good estimate. 

Plots of the LOES and the simulation time histories 
are shown in Fig. 4. The LOES matches the simulation 
data exceptionally well. The only discrepancy occurs 
between 9 and 10 sec. This discrepancy is believed 
to be the result of a nonlinearity around the zero stick 
position, possibly caused by break out forces or hys- 
teresis in the stick at that point. 

Flight 

For the flight data estimation, the final results of the 
simulation estimation are used as the initial state es- 
timates. The values for process noise are those de- 
termined from the simulation. The initial variance of 
the state is that used in the simulation estimation. The 
flight data is then processed to arrive at the final esti- 
mate, which is compared to the simulation final esti- 
mate. The output of the flight data estimation is pre- 
sented in Table 2. 

The same correlations among parameters that are 
seen on the simulation, S with Mo and 6o , and Ms with 
M v occur here. As before, this result is not desirable, 
but it is expected and is not thought to be a problem in 
the estimation. 

Plots of the LOES and the flight data are presented 
in Fig. 5 and demonstrate excellent agreement with the 
flight data. However, for reasons which are not clear, 
the filter underestimates the final peak. 

Figure 6 shows the innovations. The bounds of the 
innovations are obtained by the equation 

bound = ±2[ HPH 1 + R]t ( 23) 

The structure, especially in the pitch rate data, indi- 
cates correlation over time because a low-order model 
was used to describe higher-order dynamics, resulting 
in process noise. This correlation also suggests that 
information is still available in the signal, or that the 
noise is not Gaussian and white. Additional timing of 
the gains may reduce the correlation. 

Figure 7 shows the time histories for the state esti- 
mates. In almost all cases, the states setde by approx- 
imately 9 sec. These results are encouraging. If the 
time histories for the states representing the parame- 
ter estimates vary wildly at the end of the estimation. 


then little confidence could be held in the estimation 
for that state. Large variations would also be reflected 
in a large final variance for that state. 

Simulation and Flight Comparisons 

Figure 8 is a comparison of the final state estimates, 
with a ±2 a error bound (based on the final P), of the 
simulation and flight estimations. The frequency and 
damping ratio are computed from the state estimates; 
their bounds were computed by a simple linearization. 
This is only an estimate of the bounds and is thought 
to be sufficient because of the low correlation of the 
parameters. The results were expected to agree better 
than they do. 

The damping ratio is estimated to be higher from 
flight than from the simulation. This increased damp- 
ing ratio has also been observed by the pilots and in 
postflight frequency response analysis. This increase 
in the damping ratio is also observed in the near-real- 
time open-loop frequency response analysis as an in- 
crease in the phase margin. 

To further verify the results of the algorithm, the 
damping ratio and pure time delay estimates are com- 
pared to the results from the postflight frequency re- 
sponse analysis. The EKF estimated a time delay be- 
tween 140 msec and 146 msec; the frequency response 
analysis estimated a time delay of 120 msec. For the 
damping ratio, the EKF estimated 0.84 to 0.92, the fre- 
quency response analysis estimated 1.05. The agree- 
ment is good, considering frequency and damping ratio 
tend to tradeoff without affecting performance. 

Thus, the differences between the flight data esti- 
mates and the simulation data estimates are caused 
by a discrepancy between the simulation and airplane. 
Some discrepancies between the simulation and air- 
craft performance are expected. This demonstrates the 
EKF’s ability to identify such discrepancies. 

Robustness 

The data was also analyzed in more detail to deter- 
mine the robustness of the algorithm. In one instance, 
the window over which the estimation was performed 
was shifted 1.5 sec. This shift only has a slight effect 
on the final state estimates. 

In another check for robustness, the filter was reini- 
tialized with the final state estimate from the previ- 
ous run and the process was repeated in an iterative 
fashion until convergence. Although the results dif- 
fered slightly, they showed that the initial estimation 
would be sufficient. This is where the Kalman filter 
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possesses an advantage over other identification tech- 
niques, which require multiple processing of the data 
to obtain a final estimate. The results of the above two 
investigations are presented in Table 3. 

The most dramatic evidence of the robustness of the 
algorithm, however, is its reaction to a data spike in the 
stick input signal at approximately 15 sec (Fig. 3(a)). 
This channel is the most sensitive to data spikes since 
the model in the filter algorithm will be driven by this 
input, producing outputs inconsistent with the aircraft 
outputs. As a result large differences between mea- 
sured and predicted values can occur altering the state 
estimates unfavorably. 

The flight data was estimated in the presence of this 
spike with the process noise weighting matrix equal to 
zero. The effects of this spike can be seen most clearly 
in the plots for the filter state (Fig. 9), at approximately 
10.6 sec (the difference in time is due to the window 
over which the flight data is identified). The LOES 
model and flight data (Fig. 10) still agree well espe- 
cially considering that the process noise was zero. The 
ability of the algorithm to survive a data spike of such 
magnitude on the most critical channel indicates that 
the algorithm will perform well during flight test. 

Since the technique is in the developmental stage 
all the analysis so far has been postflight analysis. The 
goal of this technique is implementation during flight 
test. While proving the concept, however, execution 
speed was not considered. Once the technique was 
proven sound and acceptable, computation time was 
considered. By then, there was time for only a cur- 
sory investigation. The resulting execution time was 
approximately 2.5 min (central processing unit (cpu) 
time) to process 10 sec of flight data. Major sections 
of the code may be written more efficiently, reducing 
the execution time. When the program is moved to 
the computer used during flight test, faster computa- 
tion times may occur because of the increased capabil- 
ities of this system. 

Concluding Remarks 

A third-order model based on the short-period ap- 
proximation is an acceptable model to describe the 
pitch axis transient response of highly augmented and 
very large order aircraft dynamics. This equivalent 


system model can be estimated by using an extended 
Kalman filter. Results suggest such an estimation 
could be conducted in near real time during flight test. 

The results of this technique have been verified 
against the existing near-real-time open-loop fre- 
quency response analysis and postflight frequency do- 
main handling qualities analysis. 

The ability to determine equivalent time delay, un- 
damped natural frequency, and damping ratio in near 
real time provides the flight test organization with 
quick and accurate information about the aircraft, con- 
tributing to safe, quick, and efficient flight testing. 
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Table 1. Output for simulation identification. 
(Undamped natural frequency = 4.2596; damping ratio = 0.5753) 
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Tabic 2. Output for flight data identification. 
(Undamped natural frequency = 4.6965; damping ratio = 0.8823) 
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Wing flap 
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Fig. 1 The X-29A airplane. 
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c 



Pilot command 

^ P Stick signal to flight control system 

^ p Stick signal used to drive extended 

Kalman filter aircraft model 

z m Measured angle of attack, pitch 
rate, and normal acceleration 


Z c Corrected angle of attack, pitch 
rate, and normal acceleration 

z p Predicted angle of attack, pitch 
rate, and normal acceleration 

Xp Predicted state vector 

*c Corrected state vector 
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Fig. 2 Extended Kalman filter algorithm. 
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(d) Normal acceleration. 
Fig. 3 Concluded. 




(b) Pitch rate. 


Fig. 4 Low-order equivalent system response compared with simulation data (perturbation quantities) as a function 
of lime. 
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(c) Normal acceleration. 
Fig. 4 Concluded. 



(a) Angle of attack. 

system response compared with flight data (perturbation quantities) as a function 


Fig. 5 Low-order equivalent 
of lime. 
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(b) Pitch rate. 



(c) Normal acceleration. 
Fig. 5 Concluded. 
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Fig. 6 Concluded. 
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(a) Angle of attack. 

Fig. 7 Estimated state lime histories for flight data identification. 
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Time, sec 

(c) Psucdo control surface deflection. 
Fig. 7 Continued. 
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(d) La. 



(c) Mat- 

Fig. 7 Continued. 


23 


.015 



Time, sec 

(0 L s . 


900051 



(g) Mb. 

Fig. 7 Continued. 
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Fig. 7 Continued. 
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0 ) M 0 . 



(m) 8q. 


Fig. 7 Concluded. 
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Fig. 8 Simulation results compared with flight data. 
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(b) M a . 

Fig. 9 Eslimatcd stale time histories for input with data spike. 





.025 



-. 015 *- 

0 


1 


2 


3 


4 5 6 7 8 

Time, sec 

(c) L s . 


J i i l 

9 10 11 12 

900062 



(d) M S . 

Fig. 9 Continued. 
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(h) Lq. 

Fig. 9 Continued. 
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— q c , low-order equivalent system estimate 
■ — q f , flight data 
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(b) Pitch rate. 

Fig. 10 Low-order equivalent system response compared with flight data (perturbation quantities) for input with 
data spike as a function of time. 
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